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Abstract 

\ We investigate how the data from various future neutrino oscillation experiments will constrain 
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the physics parameters for a three active neutrino mixing model. The investigations properly ac- 
count for the degeneracies and ambiguities associated with the phenomenology as well as estimates 
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of experimental measurement errors. Combinations of various reactor measurements with the ex- 



pected J-PARC (T2K) and NuMI offaxis (Nova) data, both with and without the increased flux 
associated with proton driver upgrades, are considered. The studies show how combinations of 



reactor and offaxis data can resolve degeneracies (e.g. the #23 degeneracy) and give more precise 
information on the oscillation parameters. A primary purpose of this investigation is to establish 
the parameter space regions where CP violation can be discovered and where the mass hierarchy 
can be determined. We find that, even with augmented flux from proton drivers, such measure- 
ments demand that sin 2 2#i3 be fairly large and in the range where it is measurable by reactor 



experiments. 
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The worldwide program to understand neutrino oscillations and determine the mixing 
parameters, CP violating effects, and mass hierarchy will require a broad combination of 
measurements. Progress in the past associated with solving the solar and atmospheric neu- 
trino puzzles took a full suite of experiments to isolate and understand the phenomenology. 
Each additional measurement helped define the direction of future studies. One would expect 
a similar chain for the current goals where the program grows as information is obtained. 
This study attempts to see how various present proposals for next generation experiments 
(including two detector reactor and accelerator based long-baseline experiments) compare 
to and complement each other. A particular emphasis is on combining experiments to give 
improved physics parameter determination. As in the past, the best constraints on the 
phenomenology come from combining data from different processes and setups. 

I. PROCEDURE FOR THIS STUDY 

For a three active neutrino scenario, neutrino oscillations are described by six physics 
parameters: #i 3 , #12, #23, Am 2 2 , Am 2 ^, and the CP violation phase, 8. In addition, a full 
description requires knowing the hierarchy of mass state 3 relative to 1 and 2, i.e. the sign 
of Am 2 ^. (See Ref. jj], Q| for a description of three neutrino oscillation phenomenology and 
current results of global fits.) 



The oscil 
is given by 



ation probability up to second order for reactor and long-baseline measurements 



Preactor - sin 2 29 13 sin 2 A + a 2 A 2 cos 4 13 sin 2 20 12 , (1) 



Piong-baseiine — sin 2 26 13 sin 2 #23 sin 2 A 

=1= a sin 26 13 sin 5c p cos #i 3 sin 2# 12 sin 2# 2 3 sin 3 A 

+ a sin 26 13 cos<5cp cos 6\ 3 sin 26^2 sin 2# 23 cos A sin 2 A 

+ a 2 cos 2 #23 sin 2 26 12 sin 2 A (2) 

with a = Ama/Am^ and A = Am| 1 L/(4£' I/ ). 

In the investigations shown here, the full formulae for the oscillation probability have 
been used as incorporated in a computer program developed by S. Parke [3j. The higher 



order corrections for the reactor probability are quite small for the distances on the scale of 
the reactor experiments considered here, and so, approximately, a measurement of P re actor 
directly constrains the mixing parameter # 13 . On the other hand, the full expression for the 
long-baseline probability introduces many degeneracies and correlations between the physics 
parameters 823 and Sep, plus the mass hierarchy through matter effects even before experi- 
mental uncertainties are taken into account [ll J]- Therefore, a measurement of Piong -baseline 
corresponds to sizable regions in the physics parameter space. 

Of the six parameters, we assumed for this study that #12, #23) Am 2 2 , an d Am 2 3 are known 
to the precision expected from either the current program (Super-K, Minos and CNGS) or 
the future program (Nova and T2K), as shown in Table E] This leaves for determination 
813, 5, and the mass hierarchy, which are the subject of this study. 

The technique we will use for this study is as follows. We generate data according 
to the experimental expectations with the underlying physics assumptions and assumed 
parameters. We then explore the parameter space for one of the measurement parameters 
in terms of Ax 2 values. For a given set of parameters, a x 2 value is found by comparing the 
prediction with these parameters to the orginally generated data. The difference between 
this x 2 value and the one calculated with the original parameters is the Ax 2 value. In 
each case, we will explore the change due to the variation of only one parameter. However, 
estimating a correct sensitivity depends on the choices of all six of the parameters described 
above. We handle these inherent ambiguities by choosing the most conservative (lowest) 
Ax 2 value for the parameter being tested by marginalizing over all the other parameters. 
For example, if we are studying 6*13, for each value of 813 studied, the other parameters 
5cp, Am| 3 (including sign) and 6*23 are allowed to vary within the allowed ranges. This 
sometimes results in discontinuities, where one hierarchy becomes worse than the other in 
a certain region. The best choices are determined by minimizing the Ax 2 value using the 
Minuitjsl program. 

The other experimental inputs for the study are given in Table HH and are derived from 
estimates of the measurement sensitivities. For these sensitivities, we have taken the values 
from the given experiment's estimates without further study. Three types of two detector 
reactor experiments are considered corresponding to a small (Double-CHOOZ, Q|), medium 
(Braidwood 7], Daya Bay 0), or large (MiniBooNE size) detector reactor V e measurements. 
The sensitivities for the reactor experiments are scaled from the sin 2 2# 13 90% C.L. limits 
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Parameter Value 


Current Uncertainty (la) 


Future Uncertainty (1 a) 


sin 2 26*23 1.0 
Am| 3 (eV 2 ) 2.5 x lO" 3 

9 12 (deg) 32.31 
Am 2 2 (eV 2 ) 7.9 x 10~ 5 


0.1 (Super-K) 
0.55 x 10- 3 (Super-K) 
2.55 (global solar fit) 
0.55 x 10~ 5 (global solar fit) 


0.01 (T2K) 
0.1 x 10~ 3 (T2K) 



TABLE I: Current and future uncertainty estimates on oscillation parameters from Super- 
Kamiokande T2K and a global solar fit by Kamland }l~3 |. The studies here use the 
above values except the central value of sin 2 26*23 is assumed to be 0.95 (see text). 

at A772I3 = 2.5 x 1CU 3 eV 2 for a null oscillation scenario. In terms of integrated luminosity 
denned as detector mass [tons] x thermal reactor power [GW]x running time [years], the 
three options correspond to 300, 3000, and 20,000 ton-GW-yrs for the small, medium, and 
large reactor scenarios. The assumed 90% CL. limits on sin 2 26\3 for the three scenarios 
are then given by 0.03, 0.01, and 0.005, for small, medium, and large detector reactor 
experiments, respectively. 

Two offaxis long-baseline experiments are considered, J- PARC to Super-K |10J (T2K) and 



the NuMI offaxis proposal |12| (Nova). For these experiments, the uncertainties are scaled 
from the expected number of events given in the Nova Proposal for 10km offaxis [13] and 
a talk by Nakaya given at NOON 2003 14]. The given uncertainties include statistical 
errors associated with the background and signal for a 5 year data run but do not include 
any systematic uncertainty. In the studies presented here, results are given for various 
data running periods with the nominal beam rates corresponding to Table |H] In addition, 
some results are given for upgraded beams with five times the flux as would be appropriate 
for the T2K Phase II experiment or the Nova experiment with a new proton driver. We 
consider the different cases of running only with neutrinos, or with neutrino and antineutrino 
running with the number of years of neutrino and antineutrino running specified for each 
plot. The uncertainty on the #23 parameter can have a significant effect on the long-baseline 
measurements since the quantity that is constrained as given in Table H] is sin 2 2^23, and the 
parameter that modulates the long-baseline oscillation probability is sin 2 #23. This can lead 
to a 60% uncertainty in the oscillation probability for sin 2 26*23 = 0.95. 

For the studies presented in this report, the uncertainties due to the variations of #23, 
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Experiment 


Basis of 
Estimate 


Osc. Prob. and a for sin 2 2# 13 = 


0.02 


0.05 


0.10 


Reactor (E v = 3.6 MeV) 
Small 1.05 km 
Medium 1.8 km 
Large 1.8 km 


sin 2 26% mit 
0.03@90%CL 
0.01@90%CL 
0.005@90%CL 


0.013 ±0.012 
0.022 ± 0.006 
0.022 ± 0.003 


0.033 ± 0.012 
0.052 ± 0.006 
0.052 ± 0.003 


0.064 ±0.012 
0.102 ±0.006 
0.102 ±0.003 


T2K (E v = 650 MeV) 
(L) = 295 km 

V 

V 


Number of events in 5 yrs: 

105.0 signal / 17.8 bkgnd 
30.8 signal / 10.2 bkgnd 


0.010 ± 0.003 
0.009 ± 0.007 


0.023 ± 0.004 
0.020 ± 0.008 


0.044 ± 0.006 
0.038 ±0.010 


Nova [E v = 2.1 GeV) 
(L) = 810 km 

V 

V 


Number of events in 5 yrs: 

227.4 signal / 39.0 bkgnd 
109.0 signal / 18.5 bkgnd 


0.010 ± 0.002 
0.008 ± 0.003 


0.024 ± 0.003 
0.017 ±0.005 


0.045 ± 0.004 
0.032 ± 0.006 



TABLE II: Estimates of the experimental uncertainties for future oscillation experiments. The 
last three columns indicate the oscillation probability and respective error for three different values 
of sin 2 2013. For the long-baseline experiments, the number of events and sensitivity is given for 
5 years of neutrino mode and 5 years of antineutrino mode separately. The given uncertainties 
include statistical errors associated with the background and signal for a 5 year data run, but 
do not include systematic uncertainty. Am 2 3 = 2.5 x 10~ 3 eV 2 for all estimates and additionally, 
sin 2 26*13 = 0.1, 5cp = for the long baseline experiment's event rates. 

Ara| 3 , and the mass hierarchy are included. In the 17 data, there is a 5% to 10% contamina- 
tion of v induced events in the sample. For the studies presented here, we use 5% for Nova, 
and 10% for T2K, as estimated by the experiment. To summarize, the following list gives 
the parameters used for the various fit results and the values and typical uncertainties for 
the parameters that are varied: 

• Am| 3 ( ~ Am 2 3 ) = 2-5 x 10~ 3 eV 2 and allowed to vary within its future expected un- 
certainty (typical o — 0.1 x 10~ 3 eV 2 ). We include the hierarchy ambiguity associated 
with Am 2 3 being < or > 0. 
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• sin 2 26*23 = 0.95 and allowed to vary within its future expected uncertainty (typical 
o = 0.01). We include the ambiguity associated with #23 being < 45° or > 45°. 

• Am 2 2 held fixed at 7.9 x 10~ 5 eV 2 . 

• 9 12 held fixed at 32.31°. 

• Sep allowed to vary between 0° to 360°. 

• sin 2 26*13 allowed to vary between 0.0 and 0.3. 

II. 90% C.L. ALLOWED RANGE OF 13 

The #13 mixing angle is as yet undetermined in the current neutrino oscillation phe- 
nomenology. It is known to be smaller than the other mixing angles. The size of #13 is an 
important ingredient in constraining models of neutrino masses and mixing, such as attempts 
to relate the quark and lepton mixings. The size of this mixing angle also has important 
implications for long-baseline v e appearance measurements because it scales the size of the 
oscillation probability. 

As an indication of how well a given measurement can constrain the value of #13 , Figure ^ 
(left: T2K, right: Nova) shows the 90% C.L. allowed regions associated with measurements 
of a null oscillation scenario where the true value of sin 2 26*13 is equal to zero. The grey region 
(white curve) is the 90% C.L. allowed region for the two long-baseline experiments for a three 
year neutrino only run with the nominal (x5) beam rate. Combining the long-baseline and 
medium reactor measurement gives the improved black region. 

If 6*i3 is large enough, then positive signals will be observed by the experiments. Under 
these circumstances, the goal would be to make the best determination of the mixing param- 
eter. Figure El shows the 90% C.L. regions that would be obtained for an underlying scenario 
where sin 2 26*13 = 0.05. As shown, a long-baseline only measurement will not determine the 
mixing parameter 6*13 very well with an allowed region that spans from 0.025 to over 0.11. 
On the other hand, a reactor experiment with at least the medium scale sensitivity measures 
sin 2 26*13 to about 10% or Q n to ±0.4°. 

As seen from the figures, the reactor measurements are very efficient at constraining 
the value of #13 . Even a small reactor experiment can probe for an early indication if the 
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FIG. 1: 90% C.L. upper limit regions for various oscillation measurements for an underlying null 
oscillation scenario where sin 26»i3 = (sin 2 26*23 = 0.95 ± 0.01, Am 2 = 2.5 ± 0.1 x 10~ 3 eV 2 and 
&CP = 0°). The left (right) plot is for the T2K (Nova) long-baseline experiment. The grey region 
is the 90% C.L. allowed region for the long-baseline experiments for a three year neutrino only run 
with nominal beam rate. The white line is the limit of the 90% C.L. allowed region for a three year, 
neutrino only run at 5x the nominal beam rate. The black region gives the combination of three 
year long-baseline runs with a medium reactor measurement. The vertical dashed line indicates 
the 90% CL upper limit for a medium reactor experiment alone. 

value is sizable. The large reactor experiment has sensitivity comparable to planned long- 
baseline experiments and the medium scale experiment can measure values in the range 
for sin 2 2^x3 > 0.02 at the 10% to 30% level. As will be seen in later plots, studies of CP 
violation and matter effects over the next decade are only possible if sin 2 2#i3 is significantly 
larger than about 0.01. A small or medium scale reactor experiment can establish if these 
studies will be possible and, if they are, add additional information for determining the 
mixing parameters. 
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FIG. 2: 90% C.L. regions for underlying oscillation parameters of sin 2 2^13 = 0.05, sin 2 2^23 = 
0.95 ±0.01, Am 2 = 2.5 ±0.1 x 10~ 3 eV 2 and 5 C p = 0°. The grey regions are for the T2K (left plot) 
or Nova (right plot) experiments for three years of neutrino running. The black regions are the 
90% C.L. allowed regions for a combined medium reactor plus long-baseline analysis. The dashed 
lines indicate how the combined measurement would degrade with the small reactor sensitivity. 

III. RESOLUTION OF THE 6 23 DEGENERACY 

The mixing angle 823 is an important parameter in developing an understanding of the 
mixing matrix and for proceeding with a determination of #13. In many theoretical models, 
623 is not expected to be 45° and the difference from this value, both in sign and magnitude, 
may lead to a deeper understanding of the mixing. 

Information on the value of 823 has been obtained from disappearance measurements 
in the atmospheric Am 2 region, such as the Super-K and K2K experiments. These experi- 
ments restrict the allowed region of sin 2 28 2 3 . Unfortunately, a single value of sin 2 2^23 = a 
corresponds to two possible solutions for #23, \ sin -1 (yfa) or | — , sin -1 {\fa). The current 
Super-K measurement of sin 2 2# 23 = 1.00 ± 0.1 corresponds to values of #23 — 45° ± 9.22°. 
For the determination of #13 using a long-baseline — > v e appearance measurement, this 
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ambiguity presents a problem since the oscillation probability is proportional to sin 2 #23, as 
shown in Equation ^ The present Super-K measurement would correspond to a change in 
the T2K or Nova oscillation probability of about 63%, for a change in #23 from 35.78° to 
54.22°. 

This ambiguity in the determination of #23 is difficult for a long-baseline — > v e ap- 
pearance measurement to resolve, but can be well addressed with a combination of reactor 
and long-baseline measurements (see Ref. 

Figure H3 shows examples of different 
combinations of long-baseline results with (black regions) or without (grey regions) the in- 
clusion of a medium scale reactor measurement. (The dashed curve is for the inclusion of 
a small, Double-CHOOZ type measurement.) For this analysis, sin 2 2# 2 3 = 0.95 ± 0.05, 
sin 2 2# 13 = 0.05, and Sep = 270°. As shown, the medium scale reactor data resolves this 
degeneracy to some degree when combined with antineutrino running, however a small 
(Double-CHOOZ) reactor experiment does not at all. 

Future measurements of disappearance in the atmospheric Am 2 region by the T2K 
or Nova experiments could reduce the uncertainty in sin 2 2023 to order 0.01 10], but would 
still leave the #23 vs. | — #23 ambiguity. Figure 0] shows the different combinations of 
long-baseline results with or without the inclusion of a medium scale reactor measurement 
for sin 2 2^23 = 0.95 ± 0.01. Again, the medium scale reactor resolves the degeneracy when 
combined with T2K and Nova neutrino and antineutrino running, whereas the small reactor 
does not. 

In general, the regions of the 5 versus sin 2 2#i3 plane where this degeneracy can be resolved 
at the two standard deviation level is shown in Fig. El From the figure, it is seen that 
resolving this ambiguity puts a premium on including reactor data with good precision. 
Combined with the upgraded off-axis data, a medium (large) scale reactor experiment can 
cover most of the parameter space for sin 2#i3 values greater than 0.05 (0.025). If sin 2 2# 23 
turns out too much different from 1.0, resolving this ambiguity will become very important 
to making progress toward determining the CP phase and the mass hierarchy. 
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FIG. 3: 90% C.L. allowed regions for simulated data with underlying oscillation parameters of 
sin 2 26*13 = 0.05 , 9 23 = 38.54°, Am 2 = 2.5 ± 0.1 x 10~ 3 eV 2 and 5 C p = 270°. The analysis 
includes the restriction that sin 2 26*23 = 0.95 ± 0.05. The grey regions are for various long-baseline 
combinations of the T2K and Nova experiments for three year running periods. The black regions 
are the 90% C.L. allowed regions for a combined medium reactor plus long-baseline analysis. The 
dashed lines indicate a combined analysis with a small reactor measurement. 
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FIG. 4: 90% C.L. allowed regions for simulated data with an underlying oscillation parameters 
of sin 2 26»i 3 = 0.05, <9 23 = 38.54°, Am 2 = 2.5 ± 0.1 x 10~ 3 eV 2 and 5 C p = 270°. The analysis 
includes the restriction that sin 2 2#23 = 0.95 ± 0.01. The grey regions are for various long-baseline 
combinations of the T2K and Nova experiments for three year running periods. The black regions 
are the 90% C.L. allowed regions for a combined medium reactor plus long-baseline analysis. The 
dashed lines indicate a combined analysis with a small reactor measurement. 
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0.1 0.1 

sin 2 20 13 sin 2 20 13 



FIG. 5: Regions in the Sep — sin 2 29±3 plane for which the ambiguity for #23 > 45° or < 45° is 
resolved at the two standard deviation level for sin 2 2#23 = 0.95 db 0.01. (The allowed regions are 
to the right of the curves.) The left plot is the result for a combined data set using the nominal 
rate T2K and Nova data with medium (light grey) or small (black) scale reactor data. The right 
plot shows similar results for the enhanced rate (x5) T2K and Nova data combined with a small 
(black), medium (light grey), and large (grey) scale reactor data. 

IV. CONSTRAINING THE CP VIOLATION PARAMETERS AT 90% C.L. OR 
DISCOVERY AT 2a OR 3a 



One of the important goals of an oscillation physics program is to determine if CP violat- 
ing effects are present in the lepton sector, as probed through the neutrino mixing matrix. In 
contrast to the reactor disappearance probability, the oscillation probabilities for the long- 
baseline experiments are affected by the value of the CP violation phase 5cp- Due to these 
differences, combinations of long-baseline neutrino, antineutrino, and reactor measurements 
can be used to isolate these CP violating effects and place constraints on Sep (see Ref. jbsj]). 
The size of these effects is scaled by the value of sin 2 2#i3 which is therefore an important 
parameter for setting the sensitivity to CP violation. The analysis also needs to include the 
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uncertainties associated with the other parameters, especially the mass hierarchy. 

Figure El gives the — > v e appearance oscillation probability as a function of 5c p for the 
various combinations of beam type and mass hierarchy for sin 2 2#i 3 = 0.05. As seen from 
the figure, there is a dramatic difference between the T2K and Nova experiments due to 
matter effects. The goal of a combined oscillation analysis would be to use these differences 
in the oscillation probabilities to constrain the CP violation phase and the mass hierarchy. 
For these analyses, reactor measurements provide an unambiguous constraint on the value 
of sin 2 20i3. 

From Figure |B1 it can also be seen that a measurement of the appearance probability for 
neutrino running alone could give information on 5cp if the value of sin 2 2# 13 was known with 
sufficient accuracy from a reactor oscillation measurement. For high values of sin 2 2#i 3 near 
the current limit, Figure [7| shows how 5c p can be constrained with neutrino only running 
for Nova, T2K, and the combination of the two with and without a medium reactor. For 
the true value of Sep = 90°, T2K and the reactor narrows down the allowed range in 5cp- 

Figure El shows how the combinations of various measurements, including antineutrino 
measurements, can be used to constrain the allowed CP violation phase. The results are 
for a scenario with sin 2 2# 13 = 0.05 and the optimum phase point 5cp = 270° (where the 
difference between the neutrino and antineutrino oscillation probability is the largest for the 
normal hierarchy). In the upper left plot, a T2K z/-only (3 years) measurement is displayed 
first without any reactor measurement (grey region), then combined with a medium scale 
reactor measurement (black region). The upper right plot then shows what happens when 
both neutrinos (3 years) and antineutrinos (3 years) are used with and without the reactor 
measurement. Finally, the lower left plot shows the combination of T2K and Nova with 
and without the reactor result. The medium reactor measurement, in all cases, significantly 
reduces the uncertainty on 13 . 

As a measure of how well the CP phase can be constrained in general, Figure El gives the 
discovery regions in the 5c p — sin 2 2#i 3 plane for which a null CP violation solution (5c p = 
or 7r) is ruled out by at least three standard deviations. The black regions use long-baseline 
data only and the grey regions include data from a medium scale reactor experiment. Plot 
a) is for Nova only, b) is for T2K only and c) includes data from both T2K and Nova. (The 
results for nominal beam rates are not shown since they do not provide any restrictions 
at the three standard deviation level.) From the plots, it is seen that the combination of 
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FIG. 6: Oscillation probability for Un — > v e appearance vs. Sep for the T2K (top) and Nova 
(bottom) experimental setups with Am 2 = 2.5 x 10~ 3 eV 2 , sin 2 2^13 = 0.05 and sin 2 2023 = 0.95. 
The four curves correspond to pure neutrino (solid: normal hierarchy, dashed: inverted hierarchy) 
or antineutrino (black circle: normal hierarchy, white circle: inverted hierarchy) beams. 

T2K and Nova with increased intensity can start to probe the CP violation phase space if 
sin 2 29 13 > 0.02. (The narrow region in the lower range of 5 is due to the ambiguity between 
normal and inverted hierarchies.) The reactor measurements show how viable a CP violation 
measurement will be with the various combinations of the long-baseline setups. 
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FIG. 7: 90% C.L. regions for underlying oscillation parameters of sin 2 2^13 = 0.1, sin 2 2023 = 
0.95 ± 0.01, Am 2 = 2.5 ± 0.1 x 10~ 3 eV 2 and 5 C p = 90°. The grey regions are for the T2K, 
Nova or Nova+T2K experiments for five years of neutrino running. The black regions are the 90% 
C.L. allowed regions for a combined medium reactor plus long-baseline analysis. The shelf-like 
discontinuity around 5c p = 0° and 360° are due to the ability to distinguish the mass hierarchy at 
those values of 5. 
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FIG. 8: 90% C.L. regions for various combinations of oscillation results for sin 2 2#i3 = 0.05 and 
S CP = 270° (Am 2 = 2.5 ± 0.1 x 10~ 3 eV 2 , and sin 2 26> 23 = 0.95 ± 0.01). Upper Left: 3 year z^only 
T2K data with (black) and without (grey) medium scale reactor results. Upper Right: T2K v + V 
with (black) and without (grey) reactor result. Lower Left: T2K v + V plus Nova v+V with (black) 
and without (grey) reactor result. 
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FIG. 9: Discovery regions in the 5cp — sin 2 2^13 plane for which a null CP (<5c*p=0,7r) violation 

solution is ruled out by at least three standard deviations. The black regions use long-baseline 

data only and the grey regions include data from a medium scale reactor experiment, a) Nova (x5 

rate with Proton Driver) z^(3yr) + I7(3yr) data; b) T2K (x5 rate) ^(3yr) + I7(3yr) ; c) T2K (x5 

rate) ^(3yr) + I/(3yr) + Nova (x5 rate) z^(3yr) + z/(3yr) data. (Nova, T2K and Nova plus T2K 

with nominal rates is not shown since those combinations are not capable of making a CP violation 

discovery at the three standard deviation level.) The vertical dashed line indicates the 90% CL 

upper limit for a medium reactor experiment alone. 
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V. DETERMINING THE MASS HIERARCHY 



For constraining the mass hierarchy, one needs to compare measurements in a region 
where the oscillation probability changes significantly for a normal versus inverted mass 
spectrum (see Figure EJ). These type of changes can be induced by matter effects as the 
neutrinos or antineutrinos propagate through material. The Nova experiment is particularly 
important here due to the long distance the neutrinos travel through the matter of the earth. 

An accurate determination of the hierarchy is possible by combining the results from 
long-baseline neutrino and antineutrino data. Here again, the ambiguity with respect to 
the value of 5cp limits the determination to regions in the sin 2 2# 13 — 5cp plane. Figure ITU1 
shows the regions for which the mass hierarchy is resolved by two standard deviations. The 
black regions use long-baseline data only and the grey regions add data from a medium scale 
reactor experiment. 

These plots show that the mass hierarchy can be determined for limited regions with 
sin 2 2#i3 > 0.05 at the nominal beam rates. With the enhanced (x5) rates, the combination 
of T2K plus Nova covers most of the Sep range for sin 2 2# 13 > 0.03. Note that in the upper 
right plot the black and grey regions lie on top of each other, that is, a reactor doesn't 
contribute for a short joint run of Nova and T2K but does for a longer one. 



VI. OTHER STUDIES 



It may be possible in the future to use a very large detector at a site with multiple reactors 
to push the sin 2 2#i3 sensitivity down even beyond the "large reactor" scenario to the level 
of 0.003 at 90% C.L. Combining such a reactor measurement with enhanced long-baseline 
results improves coverage minimally for a CP measurement and hierarchy determination as 
compared to a medium reactor, as seen in Figure ITTI 

Figures IT2l and fTTTl show sin^cp vs sin 2 2#i3 for future T2K and Nova upgrades. For 
this study, we consider a Hyper-K upgrade data sample to be equivalent to 150 years (x20 
fiducial mass for 15 years) of T2K running and a Nova Phase II with a proton driver to be 
equivalent to 50 years (x5 increased beam rate or detector volume for 10 years) of normal 
Nova running. These plots, inspired by an earlier paper of Mena and Parke's show 
clearly the regions of degeneracy. 
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Mass Hierarchy resolved to 2a regions 




FIG. 10: Regions in the 5cp — sin 2 2^13 plane for which the mass hierarchy is resolved by two 
standard deviations. The grey regions use long-baseline data only and the black regions add data 
from a medium scale reactor experiment, a) Nova ^(3yr) + F(3yr) data; b) Nova plus T2K with 
z^(3yr) + I7(3yr) data; c) Nova (x5 beam rate) with ^(3yr) + p(3yr) data; d) T2K (x5 beam rate) 
z^(3yr) + z/(3yr) + Nova (x5 beam rate) z^(3yr) + I7(3yr) data. The vertical dashed line indicates 
the 90% CL upper limit for a medium reactor experiment alone. 



19 



Very Large Reactor vs Reactor 
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FIG. 11: Regions in the 5c p — sin 2 2#i3 plane where null CP violation is ruled out by three standard 
deviations (left plot) and where the mass hierarchy is resolved by two standard deviations (right 
plot) for a very large (VL) reactor, with sensitivity of sin 2 2^13 > 0.003 at 90% C.L. Both plots 
are for Nova and T2K (x5 beam rate) with z^(3yr) + z/(3yr) data. The black region indicates 
a measurement including the medium sized reactor, and the grey region includes the very large 
reactor. The vertical dashed line indicates the 90% CL upper limit for a medium reactor experiment 
alone. 

Ambiguities associated with the mass hierarchy produce multiple solutions and ambigui- 
ties associated with 623 give extended regions along the 6^3 direction. If # 23 = 45 degrees, the 
combination of T2K and Nova correctly determines the mass hierarchy, as shown in Figure 
[T2l and gives only one allowed parameter region. On the other hand, with 623 non maximal 
{623 = 38.54°), Figure Uni shows two allowed regions due to the 6*23 ambiguity that can only 
be resolved with the addition of a reactor measurement. 
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90 CL regions for 6 2 ~=45deg 
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FIG. 12: 90% C.L. allowed regions in the sinScp — sin 2 2#i3 plane for T2K (upper left) and Nova 
(upper right) and T2K (150yrs) and Nova (50yrs) (lower left) shown with the reactor measurement 
(in grey) and without (in black) all for equal v and V running. 023 and Am 2 are allowed to vary 
within future uncertainties. 623 = 45°. T2K and Nova combined (lower left plot) correctly resolves 
the mass hierarchy. 
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90% CL regions for 0„=38deg 
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FIG. 13: 90% CL allowed regions in the sindcp — sin 2 2^13 plane for T2K (upper left) and Nova 
(upper right) and T2K (150yrs) and Nova (50yrs) (lower left) shown with the reactor measurement 
(in grey) and without (in black) all for equal v and V running. #23 and Am 2 are allowed to vary 
within future uncertainties. 623 = 38.54°. T2K and Nova combined (lower left plot) correctly 
resolves the the mass hierarchy only with the additon of a reactor experiment. 



22 



VII. CONCLUSIONS 



Over the next decade, many experiments are planned to address neutrino oscillations 
and make improved measurements of the relevant parameters of neutrino mixing. As in the 
past, combining the results from different types of experiments with different setups will be 
necessary to map out the underlying physics. In the studies described here, we have tried 
to investigate the sensitivity using the suite of currently planned or proposed experiments. 
The complementarity of a program with several long-baseline accelerator experiments and 
several reactor measurements is clear and, if parameters are favorable, will lead to significant 
progress in the understanding of neutrino mixing and masses. 

As part of this program, reactor measurements hold the promise of constraining or mea- 
suring the 6*i3 mixing parameter and helping to resolve the ambiguity in determining the 
#23 mixing parameter. The sizes of these parameters are important inputs for models of 
lepton mass and mixing that span the range from GUTs trying to relate the CKM and 
MNS matrix to extra dimension models that have neutrinos propagating in the bulk. The 
smallness of #i 3 relative to the other angles may give a hint as to what the underlying theory 
may be. Besides leading to a better understanding of neutrino mixing, these angles, 9 13 and 
#23, are two of the twenty-six parameters of the standard model and, as such, are worthy of 
high precision measurement independently of other considerations. For # 13 , a two detector 
reactor experiment unambiguously measures the size of this angle with significantly better 
precision than any other proposed experimental technique. In addition, reactor data may 
be key for resolving the 6* 2 3 degeneracy. 

Looking towards probing CP violation and the mass hierarchy in the neutrino sector, 
the field will need several high-rate, long-baseline experiments. Here, the size of #13 will be 
important for interpreting the results and for planning a viable future neutrino oscillation 
program. In addition to setting the scale for future studies, a reactor result when combined 
with long-baseline measurements may also give early constraints on CP violation and early 
indications of the mass hierarchy. In the longer term, a combination of long-baseline ex- 
periments such as T2K and Nova will start to give some information about these effects if 
sin 2 2#i 3 > 0.05 and, with some upgrades and enhanced beam rates, will give good coverage 
if sin 2 2# 13 > 0.03. If #13 turns out to be smaller than these values, one will need other 
strategies for getting to the physics. Thus, an unambiguous determination 6*13 from, for 
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example, a medium scale reactor experiment, is an important ingredient in planning the 
strategy for this program, as well as accessing the phenomenology of neutrino mixing. 
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